arXiv:1502.03429v2 [astro-ph.GA] 19 Aug 2015 


Mon. Not. R. Astron. Soc. OOP. [T HlSl (—-) 


Printed 20 August 2015 (MN lAT^X style file v2.2) 


Internal alignments of red versus blue discs in dark matter 
haloes 

Victor P. Debattista 1 *, Frank C. van den Bosch 2 , Rok Roskar 3 , Thomas Quinn 4 , 
Ben Moore 3 , David R. Cole 1 

1 Jeremiah Horrocks Institute, University of Central Lancashire, Preston, PHI 2HE, UK 

2 Astronomy Department, Yale University, PO Box 208101, New Haven, CT 06520-8101, USA 

3 Research Informatics, Scientific IT Services, ETH Zurich, Weinbergstrasse 11, CH-8092, Zurich, Switzerland 

4 Astronomy Department, University of Washington, Box 351580, Seattle, WA 98195, USA 


Draft version on 20 August 2015 


ABSTRACT 

Large surveys have shown that red galaxies are preferentially aligned with their haloes, 
while blue galaxies have a more isotropic distribution. Since haloes generally align with 
their filaments, this introduces a bias in the measurement of the cosmic shear from 
weak lensing. It is therefore vitally important to understand why this difference arises. 
We explore the stability of different disc orientations within triaxial haloes. We show 
that, in the absence of gas, the disc orientation is most stable when its spin is along 
the minor axis of the halo. Instead when gas cools on to a disc it is able to form in 
almost arbitrary orientation, including off the main planes of the halo (but avoiding an 
orientation perpendicular to the halo’s intermediate axis). Substructure helps gasless 
galaxies reach alignment with the halo faster, but has less effect on galaxies when 
gas is cooling on to the disc. Our results provide a novel and natural interpretation 
for why red, gas poor galaxies are preferentially aligned with their halo, while blue, 
star-forming, galaxies have nearly random orientations, without requiring a connection 
between galaxies’ current star formation rate and their merger history. 

Key words: Galaxy: halo - galaxies: evolution - galaxies: formation - galaxies: 
fundamental parameters - galaxies: haloes - galaxies: kinematics and dynamics 


1 INTRODUCTION 

Intrinsic alignments of galaxies, i.e., correlations of galaxy 
shapes with each other or with the large scale density 
field, are one of the most important astrophy sical contam- 
inants for weak lensing measurements (see iTrox cl & Ishak 
120151. for a comprehensive review). In particular, coher¬ 
ent alignments of the shapes of physically nearby galax¬ 
ies can strongly bias the cosmic shear measurements 
that lie at the heart of upcoming e xperiments s uch as 
the Lar ge Synoptic Survey T elescope (llvezic et al.l l2008ll . 

Euclid ( Laurejjs_et_alk 2011_the Dark Energy Survey 

llSoares-Santos fe DES Collabora tion! 120 1211 an d the Wide- 
Field Infrared Survey Telescope ( Green et al.l 120111. If ig¬ 
nored, this intrinsic alignment bias can cause large system¬ 
atic errors in, for example, the inferred dark energy eq ua- 
tion of state (e.g. iBridle fe Kind [2 0071 : iKirk et Ml |20 1 2T) . It 
is prudent, therefore, that we develop a solid understanding 
of galaxy alignments. 
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Observationally, various forms of alignment have been 
detected over a wide range of scales. On scales ranging 
from ~ 1 to 100/i -1 Mpc, numerous studies h ave de¬ 
tecte d alignments between individual clus ters (e.g. IBinggelil 


1 19821 : IPlionisI 19941: Smargon et alj 2012lh between galax- 


et al.ll200ol: B rown et all 


Okumura et al 


20021: IHevmans et al] 


and between galaxies and their 
density (or tidal) field (e.g. 


20091: ILi et al.[|2013lh 

surrou nding large -scale 
Man dcl baum et al. 20061: 


Hiratf^et^l. [l200^JPaz_et alj 2008: Faltcnbaeher^t_aL||2009[: 


Tem pel fe Libeskindl 20131 : Tempel et al.l 20131 : IZhang et al) 

3|, 120151 1. Alignments have also been detected on smaller 


201 


scales, within individual host haloes. These include, among 
others, the alignment between brightest cluster galax¬ 
ies and their p arent cluster (e.g. ICarter fe Metcalfel Il980| : 
IBinggelil 1 19821 : Istrublel Il990h , between the orientation of 
central galaxies an d the spatial distributi o n of their satel¬ 
lite galaxies (e.g. ISales fe Lambasl l2004 iBrainerdl 120051: 
Yanget_alJ 2006^ Agustsson&Brainerdl 2006 ' L Azzaroet^jJ 


20071: piltenbacher et al] 20071: Wang et al. 2008bl T 2010l : 


Nierenber g et al. I l201ll : ILi et al l 12013! ) as well as be 
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tween the orientation of satellite 
central-satellite position vector (e.g. 


ga laxies and the 
Pereira & Kuhn 


20051: [Agustsson fc Brainerdll200fil: IFaltenbacher et al.ll20oit 


Hao et al . l [201ll : [Schneider et al.ll2013l ). An important out¬ 
come from all these studies is that the strength of these var¬ 
ious alignments depends on a variety of galaxy properties. 
Of particular relevance to this paper is the fact that align¬ 
ment effects are typically much stronger for red (early-type) 
galaxies than for blue (spiral) galaxies. 

From the theoretical side, tidal torque theory predicts 
large-scale correlations between t he angular momenta and 
shapes of dark matt e r haloes (e.g. Catelan__^J)orciani 2001; 
Catelan et al.l 1200 ll : ICrittenden et al.l l200ll : IPorciani et al l 


2002), which has be en tested and confirmed with numer 


ical s i mulations (e.g. Croft_&_Metzler[ 2000|; iHeavens et al.l 


l2000l ; Ijind 120021 ; IFaltenbacher et alJ 120081 b The tidal field 
is however too weak to directly affect the orientation 
of galaxies and must do so either through formation or 


accretion processes (ICamelio fc Lombardi 120151 ). Simula¬ 


tions have also shown that the shapes and/or spins of 


sitv distribution fe.g. Aragon-Calvo et al.l 20071: Hahn et al. 

2007c 

b|; Cuesta et al. 

20081: Paz et al. 2008: Zhang et al. 

to 

0 

0 

SSL 

Wang et al. 2011 

; Codis et al. 2013; Libeskind et al. 

2012; 

Trowland et al. 20131). Such alignments originate from 


both tidal torques, and from the preferred accretion of new 
material along the directions of sheets and filaments that de¬ 
li nea te the cosmic web fe. g. Ijin g 120021 : iBailin fe Steinmet3 
120051 : IFaltenbacher et ~ahl 120051) . On small scales, inside 
the virialized regions of dark matter haloes, non-linear 
effects such as violent relaxation and impulsive encoun- 
ters are be lieved to weaken any pre-infall alignments (e.g. 
IPorciani et al ] 20Ql). Ftowever, at the same time tidal forces 
due to the host halo may introduce new alignments, similar 
to the tidal- locking mechanism that affects the Earth-Moon 
system (e.g. Ciotti fe Duttalll994l ; lUsami fe Fuhmo~tolll997l : 
iFleck fe Kuhn] 20031 ). Indeed, numerical simulations have 
shown that dark matter subhaloes are preferentially radially 
aligned, with their major axes pointing towards the centre of 
their host halo fe.g. iKuhlen et al.ll2007l : IFaltenbacher et~ahl 
120081: IPereira et al.fl2008h . 

Hence, if galaxy shapes and/or orientations are de¬ 
termined either by the angular momentum or by the 
shape/orientation of the halo in which they form, then in¬ 
trinsic galaxy alignments will naturally emerge. The vast 
majority of galaxies, including blue, spiral galaxies but also 
the majority of red, ear l y-typ e s ystem s, are supported by 
rotation (see e.g.. lEmsellem et al] [20111. and their orienta¬ 
tion is therefore governed by angular momentum^] Thus, 
understanding intrinsic alignments requires understanding 
how the angular momentum vectors of discs align with their 
host haloes and their large-scale tidal field. Collisionless sim¬ 
ulations find angular momenta of haloes ali gned w ith their 
minor axes, with only small misalignments llDubinski 1992t 


Warrenet al. 1992; Porciani et al.l[2002: Bai lin fe Steinmet d 


20051 : Faltcnbacher ct al. 2005j). However, hydrodynamical 
simulations reveal a more complicated picture. First of all, 


1 Only the most massive ellipticals are pressure supported, and 
their orientations most likely reflect the directionality of their last 
major merger 


even in the absence of cooling, the halo and gas angula r mo¬ 
menta misalign on average by ~ 30° llv an den Bosch et al.l 
120021 : IChen et al.ll2003l : ISharma fe Steinmet J 2 OO 5 I ') . Further- 

more, when cooling is included, the formation of the disc 
modifies the shape and orientatio n of the inner halo, but 
leaves the outer halo largely intact (Kazantzidis et ai1l2()0ll ; 


IBailin et alJ 120051 : iDebattista et al. 20081 . hereafter D08). 


Consequently, the disc angular momentum is well aligned 
with the inner halo (r < O.lrvir), but is only poorly aligned 
with the minor axis at larger halocentric radii; overall hy¬ 
drodynamical simulations find that discs and haloes are mis¬ 
aligned by 30°-40° on average (ICroft et al.ll2009l : iBett et~ahl 
l20ld : lHahn et al .11 201 (T ) . This is of the same order as what is 
required to reconcile the relatively weak alignment strength 
observed, with the relatively strong ali gnments predicted 


Wan 

g et al. 2008b 

2013 

Zhang et al.l 


However, what is the cause of the colour-dependence 
of intrinsic alignment strength? It has been suggested (e.g. 
Ijoachimi et al Jfeoi .3 alibi) that blue disc galaxies are aligned 
with the angular momentum vector of their host halo, while 
red, early-types have moments of inertia that are aligned 
with those of their host haloes. The latter is motivated by 
the belief that early-type galaxies are the outcome of ma¬ 
jor mergers, combined with the idea that merger remnants 
are o riented along the direction of the last major merger 
(e.g. Ivan Haarlem fe van de Wevgaertl 1 19931 1. However, as 
already alluded to above, the majority of red, early-type 
galaxies are discs: the ATLAS 3D project finds that ~ 86% 
of early-type galaxies (def ined as galaxies that lack spiral 
arms) are rapidly rotating (lEmsellem et al.ll201 if) . Hence, we 
consider it unlikely that the origin of the different orienta¬ 
tions of red versus blue galaxies is predominantly accounted 
for by mergers. 

In this paper, we examine an alternative explanation 
that accounts for the fact that (i) dark matter haloes 


are triaxial 1 

Bardeen et al. 1986t Barnes & Efstathiou 1987; 

Frenk et al. 

19881; Dubinski & Carlberg 199ll: 

Jing & Suto 

20031; Bailin & Steinmetz 20051; Allgood et al. 

2006|), giv- 


ing rise to torques between the halo and the disc, (ii) 
dark matter haloes have substructure, which can tilt discs 
dHuang fe Carlberd 1 19971 ) and cause gravitational pertur¬ 
bations, and (iii) what distinguishe s red from blue discs is 
the lack of ongoing gas accretion. In lDebattista et~akl (|2013l . 
hereafter D13), we explored the stability of disc galaxies with 
their spins along the intermediate axis of the hal o, an orien¬ 
tatio n that has been proposed f or the Milky Wa y ( Law et ahl 
120091 : lLaw fe Maiewskil l20ld : iDeg fe Widrowl l2013h . D13 
showed that such an orientation is never stable, and discs 
can never form in this orientation even when gas angular 
momentum is along the intermediate axis. D13 also showed 
that a disc can survive for a long time off one of the sym¬ 
metry planes of the halo when gas is present, and proposed 
this as a natural interpretation for the Milky Way. In this 
paper, we extend this earlier work to address the stability of 
disc orientations in general, and the role of gas in supporting 
different orientations. We show that gas cooling provides a 
dominant mechanism for maintaining discs off the symmetry 
planes of the halo, and that this offers a natural explanation 
for the different orientations of red and blue galaxies. In 
particular, we suggest that red discs lack significant gas ac- 
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cretion, and have their orientation governed by halo torques, 
while blue discs have orientations that are set by the balance 
between halo torques and angular momentum of ongoing gas 
accretion. 


2 COLLISIONLESS SIMULATIONS 

We start by considering the evolution of disc galaxy orien¬ 
tation in the absence of complications introduced by gas. 

2.1 Constructing collisionless initial conditions 

As in D08, we form triaxial haloes via the head-on merger 
of two prolate haloes, the mselves the p roduc t of head-on 
mergers of spherical haloes dMoore et alj I2004h . The merg¬ 
ers, and all subseq uent collisi onless simulations, are evolved 
with pkdgrav llStadell 1200111 . an efficient, multi-stepping, 
parallel tree code. 

In the collisionless simulations, we use halo A of D08. 
We have verified that other halo models give similar re¬ 
sults. Halo A is constructed from two consecutive mergers. 
The first head-on merger places two concentration C = 10 
spherical haloes of mass 2.3 x 10 12 Mq 800 kpc apart ap¬ 
proaching each other at 50 kms -1 . The spherical haloes are 
generated from a d istrib ution function using the method of 
[ Kazantzidis et akl |2004|) with each halo composed of two 
mass species arranged on shells. The outer shell has more 
massive particles than the inner one, increasing the effec¬ 
tive resolution in the centre. We use a softening parameter 
e = 0.1 kpc (e = 0.5 kpc) for low (high) mass particles. As 
shown in D08, a large part of the particle mass segrega¬ 
tion persists after the mergers and the inner region remains 
dominated by low-mass particles. This merger produces a 
prolate halo. Halo A results from the head-on merger of two 
copies of this prolate halo starting 400 kpc apart at rest. 
The progenitor spherical haloes each have 1 million parti¬ 
cles, equally divided between their two mass species. The 
outer halo particles are ~ 19 x more massive than the inner 
ones. Thus, in total halo A consists of 4 million particles. 
The shape of halo A was presented in D13; Table [T] lists its 
propertie fl It has c/a ~ 0.35, but only a mild triaxiality 
dFranx et al.lll99lh T = (a 2 —b 2 )/ (a 2 — c 2 ) ~ 0.9 out to 100 
kpc. 

Once we produce the triaxial halo, we insert a disc of 
particles. The disc distribution is exponential with scale- 
length Ra = 3 kpc in all cases except model LA2d, which 
had Rd = 6 kpc, and Gaussian scaleheight Zd = 0.05-Rd- 
The discs are placed at various orientations within the halo. 
We refer to these experiments by the halo axis along which 
the disc’s spin is aligned: in ’short-axis’ (S) experiments, the 
disc spin is along the short axis of the halo, while in ’long- 
axis’ (L) experiments, the disc spin is along the halo’s long 
axis. An ’intermediate-axis’ (I) experiment has the disc spin 
along the halo’s intermediate axis. Most of the intermediate- 
axis simulations were presented already in D13; here, we in¬ 
clude a new model, IA3, with a significantly more massive 

2 We use a different convention for r2oo from D08. lYalluri et ahl 
teOlCl f and lYalluri et al.l d~20 1 2f ) , who used the radius at which 
p = 200 Pcrit- Here 7*200 is the radius within which the enclosed 
mass has average density 200 pcrit- 


Halo 

N p 

(10 6 ) 

N g 

do 6 ) 

M 200 

(10 12 M 0 ) 

7-200 

(kpc) 

b/a 

c/a 

A 

3.3 

- 

6.3 

379 

0.45 

0.35 

GP 

1.4 

1.3 

1.6 

239 

0.55 

0.55 

GT 

2.8 

2.7 

3.2 

304 

0.4 

0.32 


Table 1. The haloes used in the simulations of this paper. The 
properties listed are for the halo after the last merger and before 
the discs have been grown. N p and N g are the number of dark 
matter and gas particles within 7*200 3 and M200 is the halo mass 
within the virial radius, 7*200- Density axes-ratios b/a and c/a are 
by-eye averaged over the inner 20 kpc. 


Model 

O 

5 is: 

© 

_R d 

(kpc) 

Zd 

(pc) 

SA1 

5.2 

3 

150 

IA1 

17.5 

3 

150 

IA2 

7.0 

3 

150 

IA3 

70.0 

3 

150 

LAI 

3.5 

3 

150 

LA2 

10.5 

3 

150 

LA2d 

10.5 

6 

300 


Table 2. The collisionless simulations of this paper. M d is the 
mass of the disc, R^ is the disc exponential scalelength and is 
the Gaussian scaleheight of the disc. All models are grown inside 
halo A. 

disc. The discs are comprised of 3 x 10 5 equal-mass parti¬ 
cles. Initially the disc has negligible mass, but we increase 
this adiabatically linearly over time to a mass Md- During 
this time, the halo particles are free to move and remain in 
equilibrium while the disc particles are held fixed in place. 

Once the disc reaches the target mass, we set its par¬ 
ticle kinematics appropriate for the constant disc-height Zd 
and T oomre-Q = 1.5, as described in lDebattista &; Sellwoodl 
j200db For this, we caicuiate the potentiai using a hybrid 
poiar-grid code with the disc on a cylin dricat grid and the 
haio on a spherical grid dSellwoodl 120031) . In setting up the 
disc kinematics, we azimuthally average radial and vertical 
forces. From these initial conditions, the models are then 
evolved with pkdgrav, with timesteps refined such that 
St = At/2™ < r]{t/ag ) 1 / 2 , where e is the softening and a g is 
the acceleration at a particle’s current position. We use base 
timestep At = 5 Myr in all cases except model SA1 (which 
uses At = 50 Myr). For all simulations, we use r] = 0.2 and 
an opening angle of the tree code 6 = 0.7. Table [5] lists the 
collisionless simulations discussed in this paper0 

2.2 Briggs figures 

The orientation of a disc inside a triaxial halo is necessar¬ 
ily a 3D property of a galaxy. Consequently visualizing this 
in a single 2D figure requires a projection and some loss of 
information. Because we are only interested in the orienta¬ 
tion of the disc, which we quantify via the orientation of the 

3 The names of models in this paper are different from those in 
D08. 
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disc angular momentum, it is sufficient to indicate the two 
orientation angles of the disc angular momentum. The lost 
information is the amplitude of the angular momentum vec- 
tor, but th is quantity is of no interest here. Briggs figures 
ijBriggsl [l990l l. which are 2D polar coordinate representa¬ 
tions of the direction of vectors, are ideal for this purpose. 
A Briggs figure plots the standard two spherical angular co¬ 
ordinates 9 and <j> as the radial and angle coordinates in 2D 
polar coordinates. 

We use Briggs figures to present the evolution of the 
discs in our simulations relative to the dark matter halo. 
We hold the coordinate frame with respect to which the 
angles 9 and <j> are defined fixed, so changes in 9 and <f> are a 
result of disc tilting; because our haloes do not tumble (by 
construction) the directions of the principal axes of the halo, 
which we also indicate in the Briggs figures, do not change. 
Our long, short and intermediate axes are always defined at 
the virial radius, not at small radii where the halo shape is 
changed by the baryons. As our stellar discs are not strongly 
warped, we indicate only the orientation of the net angular- 
momentum of the stars out to a radius of 15 kpc. 

A note about how we indicate the principal axes in 
Briggs figures. We have chosen symbols in a simple visual 
mnemonic to indicate the short, intermediate and long axes: 
throughout we use a triangle, square and star, respectively 
to represent them. Thus, the ordering of the axes is the same 
as the ordering of the number of sides in their corresponding 
symbols. 

2.3 Evolution of the Short-Axis Orientation 

The top-left panel of Fig. |T| shows the evolution of model 
SA1 which starts with its disc spin along the short axis of 
the halo. The disc remains co-planar in this orientation for 
10 Gyr, despite its low mass. The net potential is every¬ 
where flattened in the same direction as the disc, even when 
only the halo potential is considered. Small tilts of the disc 
therefore lead to torques which drive damped precession but 
do not lead to a runaway tilt, making this orientation stable. 

2.4 Evolution of the Intermediate-Axis 
Orientation 

D13 showed that discs grown in halo A at the intermediate- 
axis orientation tilt very rapidly. They showed that this oc¬ 
curs even when the potential surrounding the disc is flat¬ 
tened by it to the extent that its spin is along the potential’s 
short axis. D13 argued that this indicates that the halo itself 
is unstable in this orientation. They presented evidence that 
orbits of dark matter particles that cross between the inner, 
flattened halo, and the outer elongated halo are unstable, 
which they proposed is the cause of the halo instability. 

Here, we present another simulation, model IA3, of a 
disc in halo A with a mass Ma = 7 x 10 11 Mq, i.e. four 
times more massive than the disc in model IA1, the most 
massive disc in halo A presented by D13. Fig. [5] shows the 
axes-ratios of the halo potential along the principal axes x, y 
and z, with the 2 -axis perpendicular to the plane of the disc 
at t = 0. This figure shows that x^/z^ > 1 and y^/z^ > 1, 
i.e. the high-mass disc has flattened the halo so much that 
the short axis of the halo potential is vertical to the disc all 


ro 



R [kpc] 

Figure 2. Halo only equipotential axes ratios in model IA3. The 
solid line shows x.p / z,p while the dashed line shows y&/z<&. The 
2 -axis is perpendicular to the disc at t = 0. 

the way out to 40 kpc. If the disc potential is included also, 
then the global potential becomes even more flattened along 
the disc’s spin axis. Therefore, the disc, out to large radii, sits 
in a potential which never has its intermediate axis parallel 
to the disc spin axis. Naively then this orientation may be 
expected to be stable. Indeed when we evolve the system 
with the halo particles frozen in place, the disc remains in 
this orientation for 10 Gyr. However, this orientation is not 
stable when the halo particles are live, as can be seen in the 
top-right panel of Fig. |T] As did model IA1 of D13, the disc 
tilts away very rapidly from this orientation to a short-axis 
orientation. 

2.5 Evolution of the Long-Axis Orientation 

The evolution of stellar discs placed with their spin along the 
long axis of the halo depends on the shape of the halo po¬ 
tential within which the disc is immersed. In model LAI, 
the low-mass disc is quickly driven off its original plane 
and precesses rapidly around the long axis while slowly tilt¬ 
ing towards an orthogonal orientation. This precession ends 
abruptly when the disc reaches a nearly intermediate-axis 
orientation; after this the disc drops towards a short-axis ori¬ 
entation, precessing slowly about this orientation, as shown 
in bottom-left panel of Fig. |T] The more massive disc in 
model LA2 does not tilt much, only reaching 9 ~ 15° after 
10 Gyr, but it precesses about the long axis. 

An explanation for why LAI tilts while LA2 appears 
stable can be obtained from Fig. [3] which shows the profile of 
the axis-ratios of the potential along the three principal axes, 
x, y and 2 with the 2 -axis perpendicular to the disc. Close to 
the centre of the galaxy the equipotentials are flattened like 
the disc (x^/z<s> ~ y<s>/z<i> > 1), becoming almost prolate and 
orthogonal to the disc further out (xq/z® ~ y<s>/z<s> < 1). 
The entire disc in model LAI sees an orthogonal prolate 
halo potential; the top panel of Fig. [3] shows that, for the 
halo potential, & x /$ z ~ &x/&z < 1- Thus, this orientation 
is unstable and any small perturbation of the orientation 
leads to the disc being torqued out of its original plane. In 
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Figure 1. Briggs figures (see Section 12.21 for a explanation of these figures) showing the evolution of models SA1 (top-left), IA3 (top- 
right), LAI (bottom-left) and LA2d (bottom-right). Dotted circles are spaced at 20° intervals, with the outer solid circle corresponding 
to 6 = 120° in each case. The disc spin is indicated by the big filled (black) circles at t = 0 and by the small filled (red) circles for later 
times. The small numbers besides some of the disc spins indicate the time. The open (black) star, square and triangle symbols indicate 
the direction of the long, intermediate and short axes of the haloes on large scales. 


model LA2, instead a large part of the disc is surrounded by 
a halo potential which is vertically flattened (& x /&z > 1). 
Thus, small perturbations do not lead to strong torques in 
model LA2, so it lasts in this orientation longer. We test 
this interpretation in model LA2d, in which a disc with the 
same mass as in LA2 but with a larger scalelength, = 6 
kpc, is grown in halo A. The shape of the resulting potential 
is presented in Fig. [3] this is not much different from LA2. 
However, now a smaller mass fraction of the disc is covered 


by a vertically flattened halo potential. As a result, while 
model LA2 only tilted by ~ 10°, model LA2d tilts by ~ 80°, 
as shown by the Briggs figure at bottom-right of Fig. [T] 
Fig. [4] compares the evolution of the tilt in LA2 and LA2d; 
while the tilt angle of model LA2 grows slowly, model LA2d, 
which has about twice the angular momentum, nonetheless 
tilts rapidly. This shows that it is the disc material sitting 
outside the vertically flattened potential that is responsible 
for the fast tilting in the long-axis orientation. 
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°0 5 10 15 20 

R [kpc] 


Figure 3. Equipotential axes ratios in models LAI (black lines), 
LA2 (blue lines) and LA2d (red lines). The solid lines show x$/z<$> 
while the dashed lines show y$/z$, where z > x > y are the halo’s 
long, intermediate and short axes on large scales and the z-axis 
is perpendicular to the disc. The dotted horizontal lines indicate 
an axis ratio of unity; if both axes-ratios follow this line then the 
potential is spherical. The top panel shows the shape of the halo 
potential only, while the bottom one shows the shape of the full 
potential. 



t [Gyr] 

Figure 4. Evolution of the angle between the instantaneous 
stellar angular momentum and the initial one in models LA2 and 
LA2d. Both models start in a long-axis orientation in the same 
halo and have the same mass. Model LA2d has a scale-length 
twice that of LA2. 



Figure 5. Evolution of the angle between the stellar disc angular 
momentum and the halo short axis at large radius. Solid lines are 
models with discs starting in the intermediate-axis orientation 
and dashed lines are for those starting in the long-axis orientation, 
as indicated at top-right. All models are evolved in halo A. 


2.6 Synthesis of collisionless simulations 

Fig- 0 compares the tilting of all the long- and intermediate- 
axis simulations in halo A. The intermediate-axis simula¬ 
tions most rapidly reach the short-axis orientation. At low 
disc mass, the long-axis orientation is unstable with the disc 
tilting and precessing. At high disc mass, the long-axis ori¬ 
entation is still unstable, but now the tilting rate is so low 
that the disc persists close to this orientation for most of a 
Hubble time. However, the fact that this orientation is only 
quasi-stable suggests that it may be susceptible to pertur¬ 
bations from satellites, as we demonstrate below. 

Comparing the discs during the period when the disc 
mass is growing (with the disc particles fixed in place and 
the halo particles free to move), we find that the disc in the 
short-axis orientation has a potential energy ~ 1% smaller 
(more bound) than the disc in the long-axis orientation for 
the same mass, making it the most stable orientation. In a 
fixed halo potential, it is readily apparent that the disc po¬ 
tential is minimized when the disc is in a short-axis orienta¬ 
tion; this is still true even when the inner halo is flattened 
by the disc growing within it. 


3 THE IMPACT OF COOLING GAS 


Although dark matter angular momentum tends to align 
with the halo’s short axis, cosmological simulations show 
that the gas angular momentu m is fre q uently decou ¬ 
pled from that of the halo Jvanden_^osch_et_aL 2002): 
Chenet_alJ_ 20031: iBailin et all l2005l : Sharma fc Steinmetd 


20051: iRoskar et al.H2010i ~ We now explore how the deliv¬ 


ery of gas with angular momentum along different directions 
affects the evolution of disc orientation. 


3.1 Initial conditions with gas 

We perform experiments of discs forming out of gas cool¬ 
ing within both prolate and fully triaxial haloes. All stars 
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Figure 6. Cartoons showing the merger geometries that produce 
the prolate models GPS, GP45 and GPL. The (red) filled arrows 
indicate the orientation of the gas angular momentum, while the 
(black) open arrows show the merging velocities of the spherical 
haloes. The final long axis of the merger remnant is along the 
merging direction. 



Figure 7. The shape of the dark matter density for the haloes 
of models GPL, GPS and GP45. Solid, dashed and dotted lines 
show b/a , c/a and T, respectively. The three haloes have identical 
shape. 

in these si mulations form sel f -consi stently out of the cooling 
gas. As did lAumer fe White] (l2013h . in early experiments we 
found that arbitrarily inserting rotating gas coronae within 
pre-existing triaxial dark matter haloes leads to a substan¬ 
tial loss of gas angular momentum, making the formation of 
realistic disc galaxies difficult. Our approach therefore is to 
include the gas right from the start while constructing the 
prolate/triaxial dark matter haloes via mergers. We set up 
a prolate halo with an equilibrium gas distribution by the 
head-on merger of two spherical Navarro-Frenk-White dark 


matter haloes each with an embedded spherical hot gas com¬ 
ponent containing 10% of the total mass and following the 
same density distri bution. The starting haloes are set up as 
described in iRoskar et al.l d2008l l: each has a mass within 
the virial radius (?' 2 oo — 200 kpc) of 10 12 Mg. A tempera¬ 
ture gradient in each corona ensures an initial gas pressure 
equilibrium for an adiabatic equation of state. Gas veloc- 
ities are initialized to give a spin param eter of A = 0.065 
l|Bullock et al.l EoOll : iMaccio et aljf2007h . with specific an¬ 
gular momentum j oc R, where R is the cylindrical radius. 
Each halo uses 1 x 10 6 particles in each of the gas and dark 
components. Gas particles initially have masses 1.4 x 10 5 Mq 
and softening 50 pc, the latter inherited by the star parti¬ 
cles, while dark matter particles come in two mass flavors 
(10 6 Mq and 3.5 x 10 6 M@ inside and outside 200 kpc, re¬ 
spectively) and with a softening of 100 pc. Initially the two 
haloes are 500 kpc apart and moving towards each other at a 
relative velocity of 100 kms -1 , producing prolate systems. 
We generate three such prolate haloes differing by the ori¬ 
entation of the gas angular momentum relative to the halo 
long axis. Since the remnant’s long axis is along the initial 
separation vector, we only need to incline the initial spher¬ 
ical gas haloes relative to the separation vector to ensure 
a remnant with gas angular momentum tilted relative to 
the halo long axis. We chose tilt angles of 0°, 45° and 90° 
(models GPL, GP45, and GPS, respectively). Fig. [6] shows 
the merger geometries of the three prolate haloes. The final 
haloes are quite prolate, T > 0.95, and have (c/a) ~ 0.65 
(see Fig. 0. 

We produce triaxial models by merging two copies of the 
prolate system GPS. By judicious orientation of the angular 
momenta before the merger, we are able to produce systems 
with angular momentum along either the long or the inter¬ 
mediate axes, which we refer to as models GTL and GTI, 
respectively. Model GTI was presented in D13. In models 
GTI and GTL, we first rotate GPS about the long axis so the 
angular momentum vector is along the j/-axis, then rotate 
the system about the 2 -axis by ±30°. Merging these from a 
separation of 500 kpc with a relative velocity of 100 kms -1 
produces a triaxial system with long, intermediate and short 
axes corresponding to the x, y and 2 axes, respectively. The 
merger geometries for these models are shown in Fig. [5] We 
also ran an additional model with gas angular momentum 
initially in the plane spanned by the long and intermediate 
axes (model GTE). We produce this model by rotating GPS 
about the y-axis by 90° as before, then merging two copies 
of this system both rotated about the 2 -axis by —30°. This 
merger geometry is also shown in Fig. [H] By construction, 
none of the prolate and triaxial haloes have any figure rota¬ 
tion. The shape of the three triaxial haloes is identical and 
is shown in Figure 1 of D13. 

These simulatio ns are evolved with GASOLINE 
dWadslev et alj i2004l h the smoothed particle hydrody¬ 
namics (SPH) version of pkdgrav. We use a timestep of 
10 Myr with a refinement parameter r/ = 0.175. During the 
merger, and for some time thereafter, we evolve the system 
adiabatically without any star formation. Then we switch 
on gas cooling, star format ion and stellar feedback using the 
prescriptions of lStinson et al.l (l2006l h A gas particle under¬ 
goes star formation if it has number density n > 0.1cm~ 3 , 
temperature T < 15, 000 K and is part of a converging flow; 
efficiency of star formation is 0.05, i.e. 5% of gas eligible to 
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Figure 8. Cartoons showing the merger geometries that produce 
the triaxial models GTL, GTI and GTE. The (red) filled arrows 
indicate the orientation of the gas angular momentum relative to 
the haloes, while the (black) open arrows show the merging ve¬ 
locities of the spherical haloes. The coordinate axes are indicated 
in the top panel: the long axis of the merger remnant is along the 
a:-axis while the short axis is perpendicular to the sheet (2-axis). 
The intermediate axis is along the ;/-axis. 


form stars does so per dynamical time. Star particles form 
with an initial mass of 1/3 that of the gas particle, which 
at our resolution corresponds to 4.6 x 10 4 M 0 . Gas particles 
can spawn multiple star particles but once they drop below 
1/5 of their initial mass the remaining mass is distributed 
amongst the nearest neighbours, leading to a decreasing 
number of gas particles. Each star particle represents a n 
entire stellar population with a Miller-Scalo (IMiller fc Scald 
Il979l) initial mass function. The evolution of star particles 
includes feedback from Type II and Type la supernovae, 
with their energy injected into the interstellar medium 
(ISM). The effect of the supernovae explosions is modeled 
at the sub-g rid level as a blastwave propagating through the 
ISM (IStinson et al.l I 2 OO 6 I) . We also include feedback from 
asymptotic giant branch stellar winds. We use an opening 
angle of 9 = 0.7. The timestep of gas particles also satisfies 
the condition St gas = r/cov.ranth/[(l + a)c + p^max], where 
Vcourant = 0.4, h is the SPH smoothing length, a is the 
shear coefficient, which is set to 1, B = 2 is t he vi scosit y 
coefficient and ^ m ax is described in IWadslev et alj (120041) . 
The SPH kernel is defined using the 32 nearest neighbours. 
Gas cooling is calculated without taking into account the 
gas metallicity. These prescriptions have be en shown to lead 
to re alistic Milky Way-type galaxies dRoskar et al.l 120121 . 
120131) . We set t = 0 to be the time at which we switch on 
gas cooling and star formation. 

Table [3] lists all the gas+star formation simulations. It 
includes the value of A ( Bullock et al.l 1200 ll ) , the angular 
momentum parameter of the gas. This is defined as 


Model 

Halo 

N * 
(10 6 ) 

M* 

(10 10 Mq) 

A 

GPS 

GP 

3.16 

9.4 

0.058 

GP45 

GP 

3.16 

9.5 

0.055 

GPL 

GP 

3.29 

9.9 

0.047 

GTI 

GT 

6.97 

20.7 

0.021 

GTL 

GT 

7.10 

21.1 

0.020 

GTE 

GT 

6.69 

19.9 

0.029 


Table 3. The simulations with gas and star formation. IV* and 
M* are the number and the mass of star particles at the final step. 
A is the angular momentum parameter (computed using equation 
UJ at t = 0 , when we turn on gas cooling and star formation. 


Lg I \E g \ 
GM t y M g 3 


(1) 


where L g , E g and M g is the angular momentum, energy 
and mass of the gas within the virial radius, and Mt is the 
total (gas+dark matter) mass within the same radius. As 
Fig. [U] suggests, the merger geometry should result in all 
the prolate halo models having the same gas angular mo¬ 
mentum. Instead A increases as the angular momentum di¬ 
rection changes from the long to the short axis, indicating 
that angular momentum is being transported within the gas 
corona. 


3.2 Evolution of disc orientation in prolate haloes 
with gas 

We start by considering the simpler evolution of discs form¬ 
ing in prolate haloes in which the initial gas angular momen¬ 
tum is inclined to the halo long axis by 0° (model GPL), 45° 
(model GP45) and 90° (model GPS). 

After 8 Gyr, all three systems form stellar discs of com¬ 
parable mass. Model GPS forms a stellar disc with its spin 
perpendicular to the halo’s long axis and remains in this 
orientation throughout its evolution. Within 1 Gyr the disc 
spin in model GPL settles almost into alignment with the 
long axis to within 5°, where it remains. Even when we shut 
off gas cooling at 4 Gyr, the disc still remains in a long- 
axis orientation, in good agreement with the results of the 
collisionless simulations. 

The most interesting result of these simulations comes 
from model GP45. In this model, the stellar disc never set¬ 
tles to one of the main planes of the halo. Nor does the disc 
line up with its angular momentum parallel to that of the 
gas. Instead, as Fig. [9] shows, the two angular momenta are 
misaligned by ~ 65° and remain in this orientation with¬ 
out precessing throughout the simulation. Fig. |TU] shows the 
dark+stellar mass distribution at 8 Gyr, clearly showing a 
stellar disc tilted relative to the main planes of the halo. The 
halo is misaligned with the disc already at 20 kpc, which is 
well within the region where the halo of the Milky Way is 
probed by the Sagittarius tidal stream. Modelling the Sagit¬ 
tarius stream therefore requires that the tilt of the dark 
matter halo relative to the disc also be taken into account. 
Interior to this region, the halo is flattened by the growing 
stellar disc. 

How is the disc able to remain in this orientation? We 
compute the torque on the stellar disc by first computing 
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the forces on the stars using GASOLINE in the usual way. 
The net torque on the disc is then given by 

T = Mi X F; (2) 

i 

where the sum is over all star particles and m;, i - ; and F; 
are the mass, position vector (relative to the centre of the 
galaxy) and force on the \ th star particle. We compute the 
torques on the stars from the dark matter halo ( th ) and from 
the gas (t 9 ) separately. Fig. I 111 presents the torques on the 
stellar disc at t = 8 Gyr. We also show the angular momen¬ 
tum of gas within 25 kpc at T < 10 6 K (J 9 ); this gas, which 
is destined to cool on to the stellar disc, arrives tilted with 
respect to its plane. As the gas cools from large distances, it 
feels the gravitational torque from the halo and its gains an¬ 
gular momentum along this direction. Once the gas reaches 
the region dominated by the stellar disc, its spin axis gains a 
component orthogonal to the r 9 direction, due to the torque 
it feels from the disc, before settling into the disc. The right 
panel of Fig. 1101 shows that the cooling gas reaches the cen¬ 
tre in a warped disc. The stellar disc, therefore, is absorbing 
from the cooling gas angular momentum in the direction op¬ 
posite to the direction in which the halo torque is pulling it. 
As a result, equilibrium between the gravitational torques 
and the gas inflow has been established. The equilibrium is 
not driven by the regulation of the gas inflow rate since the 
mass of the final disc is very similar in this model and mod¬ 
els GPL and GPS. Rather, the equilibrium is determined by 
the disc settling into an orientation where the tilting that 
would be induced by all torques on it and the tilting that the 
accreting gas angular momentum would cause cancel each 
other out. If the angular momentum orientation of the cool¬ 
ing gas changes with time, then the equilibrium orientation 
of the disc must also evolve. We show an example of such 
evolution below. 

Fig. []2] shows the evolution of the mass and angular 
momentum of the stellar and cool gas discs (defined as all 
gas within 50 kpc at temperature T < 10 5 K). By 5 Gyr 
the cold gas fraction of the galaxy has fallen to < 5% of the 
stellar mass. The stellar disc mass and angular momentum is 
growing at a much higher rate than the gas disc is declining, 
which is possible because gas is continuously cooling on to 
the disc. Over the duration of the simulation the cool gas 
net mass and angular momentum change is a small fraction 
of that of the stars, indicating that the gas mass reaches a 
detailed balance between inflow and loss from star formation 
and feedback. The bottom panel of Fig. [12] compares the 
gravitational torques on the stars and the rate of change of 
stellar angular momentum: 

|AL| = ILa-Lxl 

At t2 t 2 -ti ' { ’ 

The rate of change of stellar angular momentum generally 
exceeds the direct gravitational torque from the dark halo 
and gas corona in spite of the strongly prolate ( b/a < 0.7 
within 50 kpc) halo. The continuous delivery of angular mo¬ 
mentum misaligned with the stellar disc’s therefore explains 
why the stellar disc is able to persist misaligned relative to 
the halo. 

When the gas can cool, the disc reaches a dynamical 
equilibrium and does not change its orientation significantly; 
without cooling gas the disc must precess about the long 



Figure 11. Briggs figure showing the gravitational torques and 
angular momenta in model GP45 at 8 Gyr. The dotted circles 
are spaced by 20° intervals with the outer solid circle at 90°. 
The system has been re-oriented, as in Fig. uni such that the 
disc’s angular momentum, indicated by the solid (green) square, 
is along the 2 -axis (the origin), while the halo long axis (indicated 
by the filled, black triangle) is in the ( x , 2 )-plane. The solid (blue) 
circle indicates the torque, 77 ,, on the stars from just the dark 
matter halo (outside 25 kpc) while the filled (red) star indicates 
the torque, r 9 , on the stars from the gas. The open (red) circle 
and triangles represent the angular momentum of cool (T < 10 5 
K) gas within the inner 25 kpc in 10 annuli, with the central 
annulus indicated by the open circle. 


axis. An example of this precession, with cooling turned off 
and evolved adiabatically after 1 Gyr (when the stellar mass 
is 15% of its final mass) is shown in the right-hand panel of 
Fig. [9] The precession is not accompanied by a significant 
tilting away from the long axis in this prolate halo. We also 
verified that the disc precesses without tilting significantly 
when gas cooling is shut off at a later time when the stellar 
disc is more massive. 


3.3 Evolution of disc orientation in triaxial haloes 
with gas 

Having understood the role of gas cooling in establishing a 
dynamical equilibrium which maintains stellar discs off the 
principal planes in prolate haloes, we now consider the more 
complicated fully triaxial models. 


3.3.1 The Long-Axis Orientation 

The top-left panel of Fig. 1131 shows the evolution of model 
GTL. The disc forms nearly perpendicular to the halo long 
axis and remains in this orientation for the 9 Gyr of its 
evolution, in agreement with the evolution of model GPL in 
the prolate halo. 
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Figure 9. Left: Briggs figure for model GP45 showing 8 Gyr of evolution. Right: Briggs figure for model GP45 when gas cooling is 
turned off after 1 Gyr. Dotted circles are spaced by 20°. The disc spin is indicated by the filled (red and black) circles. The (blue and 
black) filled triangles mark the orientation of the hot gas angular momentum between 20 kpc and 50 kpc. The spin of both the disc 
and gas at 1 Gyr (left) and 1.5 Gyr (right) are marked by the larger (black) corresponding symbols. The long axis corresponds to the 
open (black) star. The outer solid circle corresponds to 6 = 90°; since this is a prolate (axisymmetric) halo each point on this circle 
corresponds to a short axis of the halo. 




x [kpc] z [kpc] x [kpc] z [kpc] 

Figure 10. Model GP45 at 8 Gyr with the system rotated so the disc is in the (cc, y)-plane and the halo long-axis is in the (x, z)-plane. 
Left: Disc+halo mass distribution. The contribution of the stars has been increased by a factor of 10x to make the disc more apparent. 
Right: gas mass distribution at the same orientation. 
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Figure 13. Briggs figures showing the evolution of models GTL (top-left), GTE (top-right), GTI (bottom-left) and GTI with star 
formation switched off after 3 Gyr (bottom-right). Dotted circles are spaced by 20°. The disc spin is indicated by the filled (red and 
black) circles. The (blue and black) filled triangles mark the orientation of the hot gas angular momentum between 20 kpc and 50 kpc. 
The spin of both the disc and gas at 1 Gyr (first output in GTL, top-left panel), 3 Gyr (when the disc settles in model GTI, bottom-left 
panel and star formation turned off, bottom-right panel) and 2 Gyr (coherent disc formed in GTE, top-right panel) are marked by the 
larger (black) corresponding symbols. The long, intermediate and short axes of the haloes correspond to the (black) open star, square 
and triangle symbols. The outer solid circle corresponds to 8 = 120°. 

3.3.2 The Intermediate-Axis Orientation 

D13 presented the evolution of model GT0. Briefly, D13 
showed that the disc that forms in this model never set¬ 
tles into an intermediate-axis orientation, even though the 

4 D13 refer to this model as Gil. 


gas angular momentum is along the intermediate axis. This 
tilting can be seen in the bottom-left panel of Fig. [13] The 
failure to settle into an intermediate-axis orientation occurs 
despite the fact that the inner halo surrounding the disc is 
flattened vertically like the disc, so that the disc sees only a 
potential with a short axis perpendicular to it. D13 argued 
that the intermediate axis orientation, which has been sug- 
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Figure 12. The evolution of the mass (top panel) and angular 
momentum (middle panel) in model GP45. The (black) circles 
show the stellar disc while the (red) triangles show the cool (T < 
10 5 K) gas disc. The bottom panel shows the rate of change of 
total angular momentum of the stars, indicated by the (black) 
circles, and the (full) gravitational torques on the stars, indicated 
by the (red) filled stars. 


Rested for the Milky Way j Law_gLalJl200^; Law fc Maiewskil 


l20ld : iDeg fc Widrowllioid : Vera-Ciro fc Helmi 20131 ) , is un¬ 

stable for the halo and therefore very implausible. 

The bottom-right panel of Fig. ll3l shows the effect of in¬ 
hibiting gas cooling after 3 Gyr (by which point the disc has 
already attained ~ 50% of its final mass). Then the disc tilts 
rapidly towards a short-axis orientation. Compared with the 
gas cooling case (fig. 11 of D13), the stellar disc approaches 
the short-axis orientation much more rapidly when gas does 
not cool, with a tilting rate almost twice that of GTI and 
comparable to the one in models IA1 and IA2 without gas. 


3.3.3 Model with evolving gas angular momentum 
orientation 

The top-right panel of Fig. 1131 shows the evolution of model 
GTE. In this model the mergers are set up such that the 
gas net angular momentum is initially in the plane spanned 
by the long and intermediate axes. However, torques on the 
corona from the halo and angular momentum transport re¬ 
sult in an inner (20 ^ r ^ 50 kpc) gas corona with angular 
momentum precessing slowly about the halo’s short axis. Be¬ 
tween 4 and 7 Gyr this angular momentum is within ~ 20° 
of the plane spanned by the short and long axes of the halo, 
enabling the stellar disc to grow in an almost long-axis orien¬ 
tation. Indeed at 7 Gyr the disc is perpendicular to the halo’s 
long axis. Yet once the gas angular momentum evolves away 
from this plane, it drags the stellar disc off the long-axis ori¬ 
entation by more than 20°. That gas can drive a disc off an 
equilibrium orientation demonstrates that gas cooling plays 
a dominant role in determining the relative orientations of 
discs and haloes. 


4 THE EFFECT OF SATELLITE 

PERTURBATIONS 

We now explore the effect of substructure on disc orienta¬ 
tions in the absence of gas by introducing a satellite in mod¬ 
els SA1 and LA2. We use a concentration c = 12 halo of mass 
1.3 x 10 11 Mg ( i.e. ~ 2% of M 200 for halo A). We place it 
at 150 kpc on the halo’s long axis and give it a tangential 
velocity such that its orbital pericentre is at ~ 10 kpc. In 
both cases, the satellite disrupts by 2 Gyr. The top panel 
of Fig. [TT] shows the effect of this perturbation: while SA1 
quickly returns to a nearly short-axis orientation, the more 
massive LA2 slowly tilts away from its original orientation 
for the remainder of its evolution. These small perturbations 
show that the long-axis orientation is only quasi-stable even 
when the disc is massive. The bottom panel of Fig. 1 14 1 shows 
a Briggs figure of the evolution of the orientation of model 
LA2 perturbed by the satellite. After the initial perturba¬ 
tion the disc precesses about the long axis, slowly moving 
away from that orientation. This precession indicates that 
the tilting is driven by the direct torquing from the halo, 
rather than a lingering influence of the satellite. The effect 
of the satellite therefore is only to provide the initial pertur¬ 
bation that drives the disc from the local energy minimum 
in the long-axis orientation. 

To test the effect of substructure when gas is present, 
we insert the same satellite into the prolate model GPL at 
4 Gyr, at which point the stellar disc has a mass of 6.6 x 
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Figure 15. Briggs figures showing the evolution of model GPL perturbed by a subhalo with gas cooling and star formation (left) and 
evolved adiabatically (right). Dotted circles are spaced by 20°. The disc spin is indicated by the filled (red and black) circles. The (blue 
and black) filled triangles mark the orientation of the hot gas angular momentum between 20 and 50 kpc. The spin of both the disc and 
gas at 1 Gyr are marked by the larger (black) corresponding symbols. The long axis of the halo is indicated by the (black) open star. 
The outer solid circle corresponds to 9 = 60°. 


10 10 Mg. The satellite mass represents 8% of M 200 for halo 
GP. The satellite is placed at 150 kpc on the halo’s long 
axis and has a tangential velocity such that its pericentre 
is at 20 kpc; the satellite disrupts by 2 Gyr. Fig. 1151 shows 
the evolution of the system when gas is allowed to cool and 
when gas is adiabatic. Both discs tilt away from their initial 
orientation; the case with no gas cooling tilts faster and 
further than the case with gas cooling. In the latter case, 
the disc tilting follows the angular momentum of the gas 
corona, while when gas is not cooling the stellar disc tilts 
independently of it. 

Finite perturbations from substructures will therefore 
destabilize galaxies in the long-axis orientation. When gas 
cools on to the disc, it slows down the tilting but the disc 
never returns to a long-axis orientation. Infalling satellites 
also impart angular momentum to the gas corona, further 
driving the disc off the long-axis orientation. 


5 DISCUSSION 

5.1 Cosmological effects not considered 

IValluri et al.l (120131 ) showed that the orbital structure of a 
dark matter halo surrounding a Milky Way-sized galaxy that 
forms in a fully cosmological model is not much different 
from isolated triaxial haloes, including some of the ones con¬ 
sidered here. Thus from an internal, orbital perspective, our 
simulations capture the same physics as more complicated 
fully cosmological simulations. 

There are, however, effects which our simulations do 
not capture. In the ACDM universe, dark matter haloes 


(1) host a spectrum of substructure accounting for roughly 
3 — 5% of th e total mass of the main halo itself at th e 
present time llDiemand et al. 2007 Sprmgel_et_aL|_ 2008). 

(2) slowly tumble ljDubinskilll992T Bailin fe Steinmetzll2004l : 
iBrvan fe Cresdl [20071 ) (3) and gain mass. We now consider 
whether each of these differences alter our conclusions. 


5.1.1 Halo substructure 

We have shown here that satellites lead to red discs evolving 
to a short-axis orientation, which is the lowest energy state 
of the system. The role of the satellites is to provide the ini¬ 
tial finite perturbation that allows halo torques to drive the 
disc to the short-axis orientation. For instance, a 2% mass 
subhalo initially drives the disc of model LA2 from the long- 
axis orientation by 10°, after which the torque from the tri¬ 
axial halo takes over and drives the disc further away from 
the long-axis orientation. Therefore the presence of a ACDM 
spectrum of sublialoes speeds up the rate at which initial 
perturbations build up but subs equently the evolu t ion i s 
dominated by the main halo itself, feazantzidis et ahl (I 2 OO 9 I) 
show that on average 5 subhaloes with mass > 20%Ma pass 
within 20 kpc of the Galactic Centre since redshift z — 1. 
The cumulative effect of these ~ 10° kicks is to ensure that 
the disc is never stuck in the local minimum provided by the 
long-axis orientation, instead ensuring that the disc is close 
to a short-axis orientation. 
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Figure 14. Top: disc tilting in models SA1 (solid line) and LA2 
(dashed line) after the disc is perturbed by a substructure. Plotted 
is the angle between the instantaneous angular momentum vector 
and the initial one. Bottom: Briggs figure for model LA2 showing 
the evolution after it is perturbed by the substructure. Dotted 
circles are spaced by 20° with the outer solid circle corresponding 
to 6 = 60°. The disc spin is indicated by the filled (red and black) 
circles. The halo long axis is indicated by the black open star. 


5.1.2 Halo figure rotation 


Pattern rotation leads to warped planes of stable loop or¬ 
bits. The measured pattern speeds of haloes in dark mat¬ 
ter only simulations follows a log-normal distribution cen¬ 
tred on Q p = 0.148ft. km s _1 kpc -1 . At these low tumbling 


at relatively large radius |Binney 197a: Heisler et al. 

1982; 

Maenenat 19821: Lake & Norman 1983r Purisen et al. 

1983; 

Steiman-Oameron & Pu 

risen 

9841: Martinet & de Zeeuw 

198^; Habe & Ikeuchi 

19851, 

19881: Peibel et al. 201ll). 

Peibel et al. (20lH) find that as the pattern speed sets large, 


both box orbits and long-axis tube orbits are destabilized. 
Prograde short-axis tubes are also destabilised, but are re¬ 
placed by retrograde short-axis tubes. These orbital changes 


can in princip le alter the shapes o f the haloes; however, even 
the slowest of lDeibel et al.l d201ll) ’s pattern speeds is about 
as large as the fastest pattern speeds of haloes found in dark 
matter-only simulations. Thus, in real haloes loop orbits at 
small radii should still be stable circulating about the short 
and long axes. The short-axis orientation should therefore 
still be a stable orientation for stellar discs. We therefore 
conclude that the tumbling of dark matter haloes is not 
large enough to alter our conclusions. 


5.1.3 Halo growth 

Haloes are generally believed to be continuously growing. 
For instance the virial radius and concentration of haloes 
increase even at late times. However much of this is pseudo¬ 
evolution due to the changing critical density of the uni¬ 
verse, relative to which the virial rad i us and concentration 
parameters are defined. iDiemer et al.l (l2013l ) show that for 
galaxy-sized haloes, M 200 ~ 10 12 Mq at 2 = 0, most of the 
evolution since z ~l can be attr ibu ted to pseudo-evolution 
of this kind. Likewise, IPrada et al.l (|2006l ) find little evolu¬ 
tion in the outer regions of haloes in the mass range 10 11 - 
5 x 10 12 Mq. Thus, most disc galaxies will not experience 
significant halo growth. More massive haloes, which typi¬ 
cally host elliptical galaxies, will still be growing rapidly 
and these will drive evolution not considered here. 


5.1.4 Halo axes-ratios 

One important difference between our simulations and cos¬ 
mological ones is in the shapes of the haloes. The haloes in 
our simulations are generally quite elongated, with 0.32 < 
c/a < 0.55; less than 3% of haloes of mass logM 2 oo — 
12.5 in cosmological simulations have c/a < 0.32 ; while 
> 16% of haloes have c/a < 0.55 dMaccio et al.1120071 . 120081 : 
ISchneider et al.ll2012l ). Therefore, the torques from the dark 
matter haloes that the discs in the simulations experience 
are larger than what might be expected in general. This has 
the effect that it allows gasless discs to settle into a minor- 
axis orientation faster here than in nature. On the other 
hand, gas cooling is able to tilt disc galaxies even more read¬ 
ily than in our simulations. 

We conclude that while our simulations are idealized, 
they capture most of the physics relevant for the relative 
orientation of discs and haloes. The elongated shapes of our 
haloes may bias the galaxies into stronger alignment with 
their haloes, but we already find that gas dominates the 
disc alignment in the case that gas cools on to the disc, 
which makes the gas even more d omin ant when the haloes 
are rounder than here. IPeBuhr et alj (120121 ) presented cos¬ 
mological simulations in which stellar discs were placed with 
their spin along either the short or the long axes of their 
haloes. They found that discs tumble significantly, particu¬ 
larly when they start out in the long-axis orientation, but 
generally remain aligned with one of the main axes inside 
50 kpc. These results with a cosmological population of sub¬ 
haloes and realistic halo tumbling rates are consistent with 
our findings, confirming that satellites destabilize discs in a 
long-axis orientation but do not alter the orientation of discs 
in the short-axis orientation quite so much. 
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5.2 Red galaxies 


In the absence of gas cooling, the main driver of a disc’s 
orientation relative to its triaxial halo is the stability of 
different orientations. Loop orbits are stable around the 
short and long axes, but not around the intermediate 
axis (Hcilisman & Schwarzschild 


Goodman fc Schwarzschildlll98ll; 


1979; Schwarzschild! 


Wilkinson fc James! 


1979; 


19821: 


Adams et al.ll2007l : ICarpintero fc Muzziol 2012lh This has led 

to the expectation that discs spin around either the long or 
the short axis of a triaxial halo. 

In the absence of gas, the orientation in which the disc’s 
spin and halo’s short axis align is always stable. Indeed this 
orientation has the lowest potential energy. We show that a 
disc can never remain with its spin aligned with the interme¬ 
diate axis of the halo, even if the disc flattens the inner halo 
so much that the vertical direction is the halo’s shortest axis 
throughout the disc. The equilibrium of discs orthogonal to 
the long axis of the halo is unstable, but the tilting rate is 
so low for massive discs that the orientation persists for a 
Hubble time and can be considered quasi-stable. However, 
perturbations from quite low-mass and low-concentration 
satellites will drive even massive discs away from this ori¬ 
entation, whereas a disc with its spin along the halo’s short 
axis returns to that orientation soon after a satellite per¬ 
turbs it. 

Red galaxies will therefore tend to line up with their 
spin parallel to the short axis of the halo. Indeed studies us¬ 
ing the Sloan Digital Sky Survey have found that red galax¬ 
ies have a satellite di stribution which is aligned with their 
projected major axis llSales fc Lain bas 2004; Brainerd| 20051 

rd||2 


Agustsson & Brainer d l2006i : lYang et al.ll2006l : iBailin et all 

20081 1. Alignment is not however perfect, because discs do 
not settle into a short-axis orientation immediately. 

We find that, at fixed mass, a disc in the long-axis ori¬ 
entation tilts much more rapidly when it has a larger scale- 
length. This occurs because more of the disc is exposed to 
the part of the potential that is perpendicular to the disc. 
A testable prediction of our simulations therefore is that, at 
fixed mass, red galaxies with larger half-mass radius should 
be more tightly coupled to their haloes and therefore show 
a stronger alignment. 


5.3 Blue galaxies 


The continued star formation ac tivity of blue g alax ies re¬ 
quires the accretion of gas fe.g. ISancisi et al1l2008l l. Cos¬ 
mological simulations find that gas and dark matter an¬ 
gular momentum decoupling is common fe.g. IBailin et al.l 
20051: ICroft et al.l]2009l; [Bett et al.l l20ld ; lHahn et al.ll20ld; 


Roskar et al]l2ofc 1. Our simulations show that, in order for 


the stellar disc to reach a balance between the gravitational 
torques and the angular momentum it gains from the halo 
the gas spin must in general be misaligned with the disc’s, 
regardless of w hether the gas is shock heated, or arrives 
via cold flo ws feirnboim fc Dekel 120031 : iKeres et alJ 120051 : 
IPekel et al.l^OOd b' This leads to star-forming blue discs re¬ 
maining in almost arbitrary orientation. Rather than be¬ 
ing perpendicular to the gas angular momentum, the disc 
tends to align itself at an orientation where the torque from 
the halo and the gas inflow balance. The supply of gas is 
therefore a much stronger driver of disc orientation than 


are gravitational torques from the halo, even when the halo 
globally has a quite prolate-triaxial shape. When the gas 
corona angular momentum orientation is evolving, the disc 
orientation also changes to remain in equilibrium between 
the hal o torque a nd the inflo wing g as angular momentum. 

IFraternali fc Tomassettil (l2012h estimate a star forma¬ 
tion rate of ~ 3Mgyr _1 in the Milky Way. Assuming a Milky 
Way disc mass of 4-5 x 10 10 Mg, this corresponds to specific 
star formation rate (SSFR) M,/M« ~ 0.06 — 0.08 Gyr -1 . 
While GP45 has a high star formation rate, ~ 7Mgyr~ 1 , the 
corresponding SSFR is M*/M* ~ 0.08 Gyr -1 , in good agree¬ 
ment with the Milky Way. Thus, the SSFR, and presumably 
the cooling rates, in the simulations are characteristic of real 
galaxies and we can expect that galaxies like the Milky Way 
will be off the main planes of their dark matter haloes. At 
lower gas cooling rates, and lower SSFRs, gas cooling is no 
longer able to deliver sufficient angular momentum to the 
disc and the disc drops into a short-axis alignment with the 
halo. 

Even in the presence of cooling gas the intermediate- 
axis orientation continues to be avoided. Since the distribu¬ 
tion of orientations avoids one axis, this introduces a weak 
correlation between the halo long axis and disc spin. D13 
found that the stellar disc that forms when gas angular mo¬ 
mentum is along the intermediate axis is never closer to it 
than 20°. Detailed modelling, which will be presented else¬ 
where, suggests that this tendency to avoid alignment along 
the intermediate axis results in a weak anti-correlation be¬ 
tween the major axis of the disc and the projected major 
axis of its host halo. The observed, weak alignments for blue 
galaxies therefore either require a correlation between the 
gas angular momentum and the large scale structure, or a 
contribution from galaxies in which the gas inflow rate is 
insufficient to (significantly) impact the orientation of the 
disc. 

The observed difference in relative orientation of disc 
and halo between red and blue galaxies is additional evi¬ 
dence that star-forming galaxies need to continuously ac¬ 
crete gas to maintain their star formation. Without the in¬ 
flow of misaligned angular momentum, torques from triaxial 
haloes and perturbations from subhaloes would drive the 
discs of blue galaxies to an alignment similar to that of 
red galaxies. The important difference between blue and red 
galaxies is not the presence of gas and star formation but 
the cosmic infall of gas with misaligned angular momentum 
with some of this gas cooling on to the stellar disc. A galaxy 
slowly forming stars from only a pre-existing gas disc/inner 
corona (including the recycling of gas via feedback and a 
fountain) will still tend to evolve towards a short-axis orien¬ 
tation. 


5.4 The Milky Way and Andromeda 


The shape of the Milky Way’s dark matter halo has been 
modelled numerous times via the Sagittarius stream. To 
date the best-fitting models have needed a triaxial halo, but 
have required that the disc is perpendicular to the inter¬ 
mediate axis of the ha lo (|Law_et_aL 2009i; ;Law fc Maiewskil 


120101 : iDeg fc Widrowl 120131 : Vera-Ciro fe Helm? 2013h . DR3 
showed that this orientation is particularly unstable and 
very unlikely even if the inner halo is flattened. By a pro¬ 
cess of elimination, they concluded that the Milky Way must 
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be tilted with r espect to the principal planes of the halo. 
iDeg fc Widrowl d2014l 'l. using mock data of a disc tilted rel¬ 
ative to the halo, showed that modelling under the incorrect 
assumption that the disc lies in one of the main planes of the 
halo leads to a best-fitting model with the disc perpendic¬ 
ular to the intermediate axis of the halo. This may explain 
the persistent finding that the Milky Way is perpendicular 
to its halo’s intermediate axis. 

The case for a tilted disc in Andromeda is suggested by 
th e distribu tion of a fraction of the satellites around it, which 
llbata et alJ (120131 ) showed consists of a thi n disc inclined 
relativ e to the main plane of the stellar disc. iBowden et al.l 
I 2 OI 3 I') show that the satellite disc must lie perpendicular to 
the long or the short axis of the halo if the satellite disc is 
long lived and the halo is not spherical, suggesting that the 
main stellar disc is inclined relative to the halo. 

5.5 Summary 

Our main results can be summarized as follows. 

• A stellar disc without gas can persist indefinitely with 
its spin along the short axis of a triaxial halo. Even when per¬ 
turbed by a satellite, the disc quickly settles back to this ori¬ 
entation. Instead if a stellar disc starts with its spin aligned 
with the long axis of the halo, whether the disc remains in 
this orientation or not depends on the shape of the halo po¬ 
tential. At low stellar mass, the potential surrounding the 
disc is orthogonal to it and the disc tilts towards a short-axis 
orientation. On the other hand, when the disc is massive it 
flattens the total potential and it usually can persist in this 
orientation for a long time. Perturbations by satellites how¬ 
ever permanently drive a disc off the long-axis orientation. 
Finally, not even a very massive stellar disc is stable per¬ 
pendicular to the intermediate axis of the dark matter halo, 
even when the inner halo becomes flattened. Thus, the most 
natural orientation for gasless red discs is with their spins 
along the short axis of the halo. 

• A disc with gas cooling on to it settles into equilibrium 
between the angular momentum it gains from the gas and 
the torques it feels from the halo and gas corona. The disc is 
therefore not perpendicular to the axes of the halo (unless 
the gas spin is about the short or the long axis) nor is it 
perpendicular to the spin of the gas. If the orientation of 
the gas angular momentum evolves, then the orientation of 
the disc is also forced to change. 

• These simulations therefore imply that while gas poor, 
red galaxies tend to align with their spins along the short 
axis of their haloes, the infall of gas on to blue galaxies allows 
them to linger at a wide range of orientations relative to the 
halo. Thus, stacking blue haloes lead to a nearly isotropic 
distribution. The need to avoid the intermediate-axis orien¬ 
tation introduces a weak anticorrelation between the disc 
and the halo major axis. 

• Our simulations show that the relative orientations of 
discs and haloes are set by local conditions, i.e. those within 
the virial radius of the halo. 
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